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This study utilizes aircraft measurements obtained off 
the California coast to analyze the geometric, 
thermodynamic, and microphysical characteristics of 
convective updrafts and downdrafts in the stratocumulus- 
topped boundary layer (BL). The vertical structure of 
thermodynamic and microphysical properties is consistent. 
The difference in peak droplet concentration between 
updrafts, downdrafts, and the defined environment increases 
near the cloud top. We observe larger mean droplet 
diameters in downdrafts near the cloud top, as the number of 
small droplets in cloud top downdrafts appears to decrease. 
The horizontal variability seen in the cloud droplet spectra 
and microphysics properties is likely due to aerosol 
activation and growth in updrafts, cloud top entrainment, 
and penetrating cumulus. We compute the mass flux, 0%, 
using a mass flux parameterization. Several measured 
variables produce consistent results and agree with mass 


flux calculations from other stratocumulus datasets. 
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I. INTRODUCTION 


A. IMPORTANCE OF STRATOCUMULUS CLOUDS TO WORLD CLIMATE AND 
TO NAVAL OPERATIONS 


It is estimated that marine stratocumulus clouds cover 
about 25% of the world’s oceans at any time (Kogan et al. 
1995), particularly along the west coasts of continents 
under the influence of subtropical high pressure and large 
scale subsidence inversions. The reflectivity of solar 
radiation at the cloud top is much higher than that of the 
ocean surface, therefore, solar energy reaching the surface 
is significantly reduced in the presence of stratocumulus. 
Stratocumulus clouds emit long-wave radiation’ at nearly the 
same temperature as the ocean surface due to their low 
altitude and hence do not greatly affect the infrared 
radiation emitted to space. Consequently, ‘the presence of 
stratocumulus clouds result in reduction of the net 
radiation received, and thus a net cooling of the earth 
system. (Kogan et al. 1995) estimate that an increase of a 
few percent of low cloud cover would counter the greenhouse 
warming effect, while similar decreases would double the 
warming .Thus it is apparent that, from a climate-modeling 
perspective, the understanding and accurate representation 
of these clouds in climate models is important. 

The structure and evolution of the cloud-topped 
‘boundary layer (BL) are also important to the United States 
Navy, both from a tactical and meteorological standpoint. A 
cloud-topped BL adds a level of complexity and uncertainty 
in assessing and predicting BL structure, potentially 


degrading environmental forecasts and the effectiveness of 


tactical decisions. The presence of Soest ocumalis, for 
example, significantly alters the BL inversion strength, 
which plays a role in electromagnetic and electro-optical 
propagation and the tactical evaluation of elevated ducts. 
The change in inversion strength and corresponding changes 
in cloud evolution also modifies aerosol-cloud interaction. 
The resulting changes to cloud structure therefore impact 
Navy remote sensing and surveillance capability, such as in 
the satellite analysis of ship tracks. The presence of 
stratocumulus also adds difficulty to short-term. and 
mesoscale forecasting. Consequently, an improved 
understanding of the evolution of a stratocumulus-topped BL 
could significantly enhance the Navy’s tactical use of the 


environment. 


B. PHYSICAL PROCESSES IN STRATOCUMULUS-TOPPED BOUNDARY 
LAYERS 


Figure 1 illustrates the major physical processes in 
the stratocumulus-topped BL. One of the physical processes 
is long-wave radiative cooling at the cloud top. fhe 
vertical extent of this cooling is limited to a shallow 
region adjacent to the cloud top. This strong cooling 
generates cold downdrafts descending from the cloud top and 
compensating warm updrafts, resulting in positive buoyancy 
flux near the cloud top. This radiative cooling is therefore 
an important source for turbulence. Another important 
process in the cloudy BL is entrainment near the cloud top, 
where turbulence in the BL mixes air from above the 
inversion into the BL. Because the entrained air is 
typically warmer and drier than BL air, entrainment mixing 


may result in evaporation of the cloud droplets and thus 











evaporative cooling: This process, under certain 
circumstances, modifies the BL and cloud structure as 
discussed in Randall (1980), Siems et al. (1990), Wang and 
Albrecht (1994), and others. 
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Figure 1. Physical Processes in the Marine Stratocumulus- 


topped BL. 


Heating from absorption of short-wave solar energy also 
occurs in the cloud layer. In contrast to the long-wave _ 
cooling, this heating is distributed more evenly in the 


cloud layer (Fravalo et al. 1981). In addition, infrared. 


(IR) radiation also results in a small amount of warming 
near the cloud base (Fravalo et al. 1981). The net result of 
both solar and IR radiation heating is to warm the cloud 
layer relative to the sub-cloud layer, forming a stable 
layer near the cloud base. . This leads to a separation of 
the upper BL (including cloud layer) from the layer below, 
prohibiting vertical transport between the two layers. This 
phenomenon is referred to as decoupling.” When decoupling 
happens, the supply of moisture to the cloud is reduced, and 
enhanced thinning of the cloud may occur. . 

The ocean surface is always a source of water vapor to 
maintain the cloud layer. However, the buoyancy effect of 
the ocean surface on the cloud-topped BL is rather uncertain 
because of the small air-sea temperature difference. The 
following chapters will present results from different 
flights that show the variability created by this effect. In 
addition to buoyancy, vertical shear of the mean wind speed 
near the surface and the cloud top modifies the dynamics of 
the cloudy BL. 

The dynamic processes within the stratocumulus cloud 
can cause vertical and horizontal variations in the droplet 
spectra. For example, in regions of dry air entrainment and 
associated downdrafts at the cloud top, droplet 
concentrations may be less than in the updraft regions, due 
to evaporation. Many studies (Rogers and Telford 1986, 
Nicholls 1989, Khalsa 1993, and Wang and Lenschow 1995) have 
found narrow regions of lower than average liquid water 
content near the stratocumulus cloud top, attributed to 
entrainment of dry air. Gerber (1996) calculated the 
effective radius of droplets inside and outside regions of 


depleted liquid water content near the cloud top and found 


little change in the relative size distribution between the 














two regions. This finding supports the suggestion of 
Nicholls and Leighton (1986) that entrainment at the cloud 
top does not effect the shape of the spectrum, but only acts 
to dilute the droplet concentration. 

Conversely, due to radiative cooling and other 
processes at the cloud top, pockets of high supersaturation 
can exist, creating a region of higher cloud droplet 
concentration (Kogan et al. 1995). Updraft regions in 
decoupled BLs often show a bimodal distribution, especially 
near the cloud base, as air parcels descending from the 
cloud top are strongly decelerated and can be recaptured by 
updrafts near the cloud base and recycled (Kogan et al. 
1995): 


C. OBJECTIVES 


The general dynamics of marine stratocumulus have been 
extensively researched, but the dynamics and microphysics in 
relation to the convective activities in these clouds are 
relatively unexplored. One subject is the role of in-cloud 
updrafts and downdrafts in determining the horizontal 
structure of the cloud and in altering droplet distributions 
within the cloud. The vertical motions within cumulus 
clouds are roughly an order of magnitude larger than in 
stratocumulus clouds, but even though the vertical motions 
are weaker, there are apparent differences between 
microphysical and dynamical parameters in stratocumulus 
updrafts and downdrafts. Understanding the relationships 
between cloud microphysics, dynamics, and radiative 
properties is an important step in the accurate modeling of 
stratocumulus development, with cloud updrafts and 


downdrafts being the important ank between these 


properties. This understanding is dependent on the 
acquisition and analysis of data capturing these important 
quantities, but few observational studies have made these 
necessary comprehensive measurements. | 

Most numerical models of stratocumulus clouds employ a 
bulk parameterization of cloud microphysics, rather than 
explicit resolution of cloud droplet size distributions. 
The former saves a large amount of computer time; but at the 
expense of accuracy and loss of potentially important 
physical effects. Advancement in microphysics formulations 
may be critical for improving prediction of stratocumulus © 
evolution. Kogan et al. (1995) employ explicit formulation 
of microphysical processes in a large eddy simulation model 
of stratocumulus, and demonstrate that cloud microphysical 
parameters are significantly affected by cloud dynamics, and ) 
that the cloud microstructure is significantly asymmetric 
between updrafts and downdrafts. The resultant mixing due to 
the updraft/downdraft plumes causes variability in the cloud 
layer structure, especially near the cloud top and base. 

This study intends to describe the significant 
horizontal variability associated with the turbulent 
updrafts and downdrafts within the stratocumulus-topped BL. 
We will analyze the characteristics of the convective events 
and the role of updrafts and downdrafts on both the 
thermodynamic and microphysical structure of the BL using a 
comprehensive data _ set. This analysis will provide 
additional understanding of the relationships between the 
cloud microphysical structure and the dynamic structure, for 
improvement of numerical modeling of these clouds in large 
eddy simulations and higher order closure models. 


This thesis is organized into five chapters. Chapter II 


describes the experiment and the data used in the analysis. 








Chapter III provides a detailed analysis of updraft and 
downdraft events within the stratocumulus-topped BL. Chapter 
IV discusses the microphysical and thermodynamic properties 
within the updrafts and downdrafts, and discusses’ the 
validity of mass flux parameterization in the stratocumulus- 


topped BL. Chapter V contains the summary and conclusions. 
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II. OVERVIEW OF THE EXPERIMENT AND THE DATA 


A. OVERVIEW OF THE EXPERIMENT 


During June and July 1987, the First International 
Satellite Cloud and Climate Project (ISCCP) Regional 
Experiment (FIRE) Marine Stratocumulus Intensive Field 
Observations (IFO) was conducted off the coast of southern 
_ California. FIRE obtained a comprehensive set of 
measurements using 12 coordinated multi-missions with five 
separate research aircraft, surface-based measurements, and 
satellite imagery. One of the main objectives of FIRE was 
to provide observations that would test theories on the 
development, ReineensAce, and dissipation Of marine 
Sstratocumulus clouds. The knowledge gained could be used to 
improve cloud parameterizations used in large-scale models 
and to develop and test satellite cloud-retrieval schemes. 

Albrecht et al. (1988) and Cox et al. (1987) give a 
detailed description of project objectives. In addition, 
Kloesel et al. (1988) gives a summary of experiment 
operations and general meteorological conditions during the 
period. A later section will briefly describe each flight 
analyzed in this thesis. In this study, we utilize data 
obtained by the NCAR Electra aircraft in seven of the ten 
flights flown during FIRE to analyze various aspects of the 
spatial variability in-the cloudy BL. The seven flights 
chosen to study contain horizontal turbulence leg data at 


tJevels boch 1m and out-of the cloud. (flights: 2, .3; 4). 6, 7; 


8 and 9). McDowell (1999) gives the general locations of 
each of the ten Electra flights. 


B. AIRCRAFT AND INSTRUMENTATION 


The National Center for Atmospheric Research (NCAR) 
Electra aircraft made turbulence, radiation, and cloud 
microphysics measurements in the BL. Each Electra mission 

consisted of horizontal legs at different altitudes in and 
above the marine BL, and vertical soundings throughout the 
BL. The aircraft speed was held constant at 100. mM: Ss": 

Horizontal flight legs consisted of 20-minute L-shaped legs 
at levels in, above, and below the cloud layer, and one 
level near the surface. These legs captured cloud structure 
and measured the turbulent flux profiles at various levels 
in the BL. Slant path or spiral soundings at multiple time 
periods and locations provided vertical profiles of BL 
variations, inversion layer characteristics, and a 
representation of the air properties above the BL.’ These 
soundings were treated as vertical profiles assuming near 
horizontal homogeneity. A typical flight track of the 
Electra measurement is shown in Figure 2. 

Measurements utilized in this thesis work include 
cloud droplet size distribution, temperature, horizontal 
and vertical wind velocity, water vapor specific humidity, 
liquid water content, and ozone concentration. Most of the 
data was sampled at a 50 Hz sampling rate and recorded at 
20 Hz for data processing. The 20 Hz data are used for 
analysis in this study. McDowell (1999) describes the 


Electra instruments used in obtaining these data 
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Figure 2. Example from FIRE Flight 2. Left: the flight 
path of the NCAR Electra during the horizontal turbulence 
legs. Right: the time variation in aircraft altitude for 
one flight segment, showing the order of the vertical 
soundings and horizontal turbulence leg measurements. 
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except for the ones described below. 
The Particle Measuring Systems (PMS) Forward | 
Scattering Spectrometer Probe (FSSP) measures concentration 
and droplet size in the 3-45 micron range. The FSSP is an 
optical particle. «counter (OPC) that detects Single 
particles and sizes them by measuring the intensity of 
light that the particle scatters when passing through a 
light beam. The Mie scattering theory 1s used to relate the 
intensity to the particle size. The size is categorized 
into one of 15 channels, each covering a 3 micrometers 
range. This study utilizes this cloud droplet size 
distribution to generate droplet spectra at different 
“yertical levels and in updrafts and downdrafts within the 
cloud. | 
A Rosemount 858 probe, a flow angle sensor that 
accounts for aircraft attack angle and sideslip, measured 


the wind velocities used in this study. 


C. DATA SELECTION 


Horizontal flight legs at different levels within and 
below cloud were used to analyze the variation in BL and 
cloud characteristics at different altitudes within the 
cloud. In general, the horizontal turbulence legs were 
identified using aircraft altitude and heading data. A 
horizontal turbulence leg was defined for each time segment 
when the aircraft was at both a constant altitude and 
heading. The in-cloud legs were identified as time 
segments within those horizontal legs where the aircraft 


measured a significant droplet concentration. Flight log 


TZ 














entries and soundings done in the vicinity of the 
horizontal legs were used to determine the cloud boundaries 
and BL heights used to scale the altitudes of the 
horizontal measurements to be presented in the later 


chapters. 
D. SUMMARY OF FLIGHT CONDITIONS 


; The following section will give a brief description of 
the BL meteorological conditions of the seven FIRE flights 
analyzed in this study based on information from Kloesel et 
al. (1988) and lLaufersweiler and Kloesel (1991). The 
general turbulence structure in each flight will also be 
discussed. These discussions will be important in 
explaining the differences seen in the cloud dynamics and 


microphysics in the data analysis contained in later 


chapters. 


1. Meteorological Conditions 


The following is a brief summary of the synoptic and 
BL meteorological conditions for each of the FIRE flights 
analyzed in this study. Table 1 lists the flight date, 
time, average BL depth, maximum liquid water content, and 
cloud depth for each flight. 

Measurements of flight 2 were taken from late morning 
to mid-afternoon. An intensified Pacific subtropical high 
to the north and an inland thermal low created 
northwesterly winds in the project area. The winds in the 

i 


BL averaged 3 to 5 ms", relatively weak in comparison to 


the other flights. The flight began with large cloud areas 


de: 


with some breaks, but by the end of the mission, a nearly 
solid cloud deck was observed in the operation region. The 
horizontal in-cloud turbulence legs were flown in a mix of 
thin, broken, and solid cloud. Broken cumulus was also 
observed beneath the stratocumulus deck, with some cumulus 
extending upwards into the stratocumulus. Drizzle was also 


observed at various times during the flight. 


Table 1. Flight and Boundary Layer Information for the 
Seven Flights Analyzed in this Study. 


Start/ Average 


End Time 


(UTC) 


6/30/87 | 1854/2348 896 
7/3/87 | 0116/0604 890 


Tio ot 1743/2304 





*Widely varied due to complex cloud conditions. 

Measurements for flight 3 began in late afternoon. A 
relatively weak Pacific subtropical high and weak thermal 
low east of the Baja Peninsula created generally northerly 
winds (7 to 10 m s') that were much stronger than for 


flight 2. The BL in the measurement area was covered by 
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nearly solid cloud with some thin breaks, with a large 
clear area to the west. As measurements continued well into 
hours of darkness, the cloud deck became very solid. The 
bases of the stratocumulus were ragged and the tops were 
variable and undulating. No lower cumulus layer was evident 
and very little drizzle was observed. 

Measurements for flight 4 started in the late morning 
and continued into the mid-afternoon. The persistent 
Pacific subtropical high and inland thermal low created 
rather strong north-northwesterly surface winds (8 to 13 m 
s*) in the measurement area on this day. A thin, solid 
‘stratocumulus deck existed, which became more broken later 
in the flight. No precipitation or lower cumulus was 
evident. | 

Flight 6 took place off of the central California 
coast. Measurements started around local noon = and 
continued into the mid-afternoon. The Pacific subtropical 
high had moved slightly westward, weakening the pressure 
gradient over the measurement area and creating moderate 
north-northwesterly surface winds. A mix of thin, solid 
and broken stratocumulus existed. Later in the flight, the 
stratocumulus became more solid, the inversion was 
relatively higher, and the cloud bases became variable, 
with some reported cloud bases as low as 100 m. All of the 
in-cloud horizontal turbulence legs analyzed in flight 6 
were from this later segment. No precipitation or lower 
cumulus was evident during the turbulence legs. 

Measurements for flight 7 started around 1100 local 
time and continued into the mid-afternoon. The Pacific 


subtropical high remained in nearly the same position as 


iis 


the previous day, while the inland thermal low weakened. 
This created a weaker pressure gradient and thus weaker 
winds out of the north-northwest in the measurement area. 
Cloud conditions were similar to flight 2 with a mix of 
thin, broken, and solid cloud, and an extensive lower 
cumulus layer with some cumulus extending upwards into the 
stratocumulus. The lower cumulus appeared to be suppressed 
in areas where the stratocumulus was thin or broken. 
Numerous showers were observed during the turbulence legs. 

Flight 8 also started in _ the late morning and 
continued into the mid-afternoon. The subtropical high had 
moved farther north and east than in previous flights and. 
the inland thermal low strengthened producing strong north- 
westerly winds (7 to 10 m s ) -in the project area. A 
relatively solid and thick cloud deck covered about half of 
the BL depth, which became slightly more broken later in 
the flight. A lower cumulus layer also existed, which in 
some cases extended up into the stratocumulus deck. Some 
drizzle was also observed. 

Measurements for flight 9 started around 0830 local 
time and continued until around noon local time. Although 
the aber oouesil high and the inland thermal low both 
strengthened from the previous day, the pressure gradient 
and winds over the project area remained relatively 
unchanged. Flight 9 was conducted in the vicinity of San 
‘Nicholas Island within the southern California bight. All 
other flights were conducted much farther from land over 
the open ocean. This was likely a major factor in the 
dramatic differences in BL cloud conditions seen in this 


flight, which were by far the most complex of all the FIRE 


LG 














flights. At times three cloud decks existed: two 
stratocumulus decks at different vertical levels, varying 
from solid to broken, and a lower cumulus field. Showers 
and fog were also observed. 

It is seen from Table 1 that all measurements, except 
flight 3, were made during the daytime. The BL height, the 
cloud depth, and the maximum liquid water content are 


comparable in all flights. 


2. Boundary Layer Mean and Turbulence Structure 


An example of the BL vertical structure from aircraft 
soundings is shown in Figure 3. This sounding is typical of 
the stratocumulus-topped BLs analyzed in this study. Here, 
the BL, as indicated by the virtual potential temperature 
profile (Figure 3a), is approximately 800 m deep, capped by 
a sharp temperature inversion. The number concentration of 
droplets (Figure 3c) remains fairly constant throughout the 
roughly 500 m depth of the cloud deck. The average diameter 
of droplets (Figure 3d) increases with height within the 
cloud, and liquid water content (Figure 3b) increases 
smoothly with height, reaching a maximum at the cloud top. 
Note that readings of droplet size below the cloud base are 
due to instrument noise and are therefore erroneous. 

Table 2 lists values of the surface buoyancy flux, 
surface frictional velocity (u+), and the cloud layer 


convective velocity scale (Wc), which are used to 


cy 
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Figure 3. Typical Vertical Profiles in Stratocumulus- 
topped Boundary Layers. (a) virtual potential temperature 
(ayy (bo) liguid water content (g kg he (c) droplet 


concentration (cm), and (d) average droplet diameter (pm). 
The measurements are from FIRE flight 9 on 16 July 1987. 
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characterize the source of turbulence in the BL. The 
convective velocity is defined using the buoyancy flux 
profiles in the cloud mixed layer following the method in 


Deardorff (1980): 
We, = {2.5(g/0,) fw Ody, (1) 


where @, is the average virtual potential temperature in the 


layer of integration, w@’ is the buoyancy flux calculated 


from direct eddy correlation method (Stull 1988), and g is 
the gravitational aeceleration. The lower limit of the 
integration, Zo, is the altitude at which buoyancy flux is 
zero which is assumed at the base of the cloud mixed layer. 
The upper limit, zi, is the cloud-top height. These limits 
are similar to those used by Deardorff, except that his 
lower limit, -zo, is the height of the cloud mixed layer base 
for decoupled cases, and O for well-mixed cases. The 
convective velocity thus defined indicates the forcing due 
to radiative cooling and other processes at the cloud top. 


The frictional velocity is defined as (Brost et al. 1982): 


2 2 
u, =(w'y’ + y'y’ ae ee 


— ee 





where wv’ and wu’ are the momentum fluxes for the north- 
south and east-west wind components,, respectively. ux 
therefore denotes the magnitude of the wind shear near the 
surface. 

The extent of vertical turbulent mixing is 
qualitatively defined as a decoupled or a well-mixed BL in 


Table 2. The definition is based upon the comparison 
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between the lidar-measured cloud-base height and _ the 
LALting condensation level (LCL) calculated using 
measurements from the lowest flight legs. In cases where 
the lidar cloud base height was unavailable, signatures of 
layer separation from the sounding profiles are used. The 
results in Table 2 for flights 2, 3, 4, and 8 are from Wang 
(1993). The same method was used for the rest of the 
LLights. Flight 7 was classified as decoupled, since the 
LCL and cloud-base height were not in good agreement. 
Since the lidar was inoperable during flight 6, soundings 
were inspected for evidence of a transition layer. The 
vertical profiles of moisture and potential temperature 
indicate a well-mixed BL, and the LCL and cloud base 
heights estimated from liquid water and droplet 
concentration jumps in sounding data are in fairly good 
agreement. The BL structure for flight 6 will therefore be 
Classified as well-mixed. FlTont 9 illustrates the 
potential variability and complexity of near-coastal 
stratocumulus-topped BLs. The complexity is suggested by 
the BL cloud structure as described earlier. There is 
significant scatter in the TKE and turbulence flux 
throughout the BL. The low value of surface buoyancy flux 
is a result of the averaging of both positive and negative 
near-surface flux values, indicating significant horizontal 
variability in the near-surface conditions. This scatter, 
coupled with the relatively smaller amount of data 
collected during the horizontal turbulence legs, makes 


classification of the BL during flight 9 difficult. 
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Table 2 Buoyancy and Shear Forcing and Mean Turbulent 
Kinetic Energy (TKE) for each Flight. The surface buoyancy 
flux is in W m%, surface friction velocity is in m Ss 
cloud layer convective velocity is in m s', the DLOCUcCEZON 
rate of near surface buoyancy and shear-generated TKE are 


in m? s? and BL averaged TKE are in m ss”. 
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* Structure of turbulent mixing is not clearly defined. 


Following Wang et al. (1999), the upper BL associated 
with the cloud layer will be referred to as the cloud mixed 
layer, while the lower layer will be called the surface- 
based layer. 

Vertical profiles of TKE from all seven flights are 
shown in Figure 4 to aid the discussion of turbulent 
structure. The BL mean TKE is given in Table 2, and 


indicates that two of the well-mixed cases (flights 3 and 


ZA 





Figure 4. Vertical profiles of TKE for all Seven Flights. 
The horizontal axes are TKE in m’s*, the vertical axes are 
height, z, scaled by the boundary layer depth, zi (Table 
1). The dotted line shows the average height of the cloud 


base. 
22 














Figure 4, Continued. 
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4) have the largest BL mean TKE. This is particularly true 
in the nocturnal case of flight 3 which also has the 
largest buoyancy forcing from the cloud mixed layer 
compared to all the other daytime flights. In the decoupled 
cases, the TKE profiles show rather weak turbulence in the 
cloud mixed layer compared to the auptace—pased layer. 
However, the convective velocity scale (ws-), or the eloud= 
top buoyancy forcing, is not particularly small relative to 
the well-mixed cases (except for flight 3). These 
observations appear to suggest that the surface forcing is. 
crucial in maintaining turbulence in the upper BL. 
Decoupling thus not only cuts off the moisture supply to 
the cloud layer, but also cuts off the TKE transport from 
the surface-based layer to the cloud-mixed layer. Further 
in-depth comparison with more well-mixed BLs is necessary 
to consolidate this finding. 

Comparison of columns 2 and 3 in Table 2 show the 
relative importance Of the buoyancy and shear, 


respectively, in generating turbulence. The near-surface 


buoyancy-generated turbulence is calculated as rida and 


the near-surface shear-generated turbulence is calculated 
as —(wu gm —). Both are terms of the TKE budget 
IZ | _ 


equation (Stull 1988). 

‘We found substantial variations from flight to flight 
in the surface forcing. The wind shear and the turbulence 
generated are significantly stronger in flights 6 and 8 
compared to other flights, and are weakest in flight 2. 


The buoyancy-generated turbulence varies from -2.34 m’s* to 
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5.10 m’s’°. The strong variation of sea surface temperature 
in the coastal upwelling region is the main reason for this 
large variation (e.g., Rogers et al. 1998). In all flights 
except flight 2, the shear-generated turbulence is much 
larger than the buoyancy-generated turbulence, which is 
consistent with eastern boundary near-surface conditions 


being close to neutral. 


Zo 
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III. GEOMETRIC PROPERTIES OF CONVECTIVE EVENTS 


A. CONDITIONAL SAMPLING METHOD 


This study uses a conditional sampling method to identify 
the turbulent updrafts and downdrafts responsible for 
vertical transport within the mixed layer and entrainment at 
the BL top. Nicholls (1989), Lenschow and Stephens (1980), 
Khalsa and Greenhut (1985), Khalsa (1993), and Wang (1993) 
use a Similar method for both the clear and cloudy BL. 
Conditional sampling is a technique in which a time series 
is selectively sampled according to specific criteria. The 
criteria are applied to a selected indicator variable. 
Vertical velocity is commonly -chosen as ‘the indicator 
variable, in which case updrafts and downdrafts within the 
BL are identified. li water vapor is used as an event 
indicator variable, as in Lenschow and Stephens (1980), the 
selected event will be moist/dry events. The event 
indicator variable also defines the boundaries of each event 
so that the geometric characteristics of the events can be 
identified. The criteria for event selection therefore 
include both the magnitude of the indicator variable and the 
size of the event. Khalsa (1993) describe the technique of 
conditional sampling in detail. 
Wang (1993) tested different thresholds of vertical 
velocity and event width in selecting appropriate criteria 
| for the conditional sampling of the FIRE data. She found 
that the results using the thresholds from Nicholls (1989) 
were very similar to those using thresholds based on visual 
inspection of the FIRE data. This suggested that the 


Nicholls criteria captured most of the convective elements 
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significant to turbulent transport in the FIRE region. 
Based on these results, this study will use conditional 
sampling criteria similar to Nicholls (1989). 

The conditional sampling scheme used in this study is 
described briefly here. The vertical velocity data is band- 
pass filtered prior to the application of the criteria to 
eliminate variations at scales smaller than 40 m and larger 
than 3 km. This prevents selection of small spurious events 
while at the same time preventing larger events from being 
discounted due to small breaks. The boundaries of an event 
are determined by the zero-crossing point of the filtered 


vertical velocity. Based on the Nicholls’ criteria, 


threshold vertical velocity is ee ws being the cloud 


convective velocity scale, and the minimum length of the 


event is 4p Zi being the BL depth. Tables 1 and 2 show 


that the BL height from all flights are similar, while three 
of the flights had a ws between 0.3 and 0.4 ms”. For 
simplicity, the length criterion is set at a constant 40 m, 
and the vertical velocity criterion is set at a constant 0.2 
ms! for all flights. The event selected will be tagged by 
a new time series called the indicator function, where the 
updraft, downdraft and background environment are 
represented by 1, -1, and 0, respectively. Figure 95 shows 
an example of the identified events. 

This study will analyze the updraft and downdraft 
geometric characteristics and will present the results in 
the next section. Chapter IV will investigate the leg- 
averaged properties of thermodynamic and microphysical 
quantities within the updrafts and downdrafts. The objective 
is to understand. the convective activity in the 


stratocumulus-topped BL, and specifically the role of 
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Figure 5. Examples of the Conditionally Sampled Updrafts 
and Downdrafts Identified Using the Identification 
Functions of 1, -1, and O for Updrafts, Downdrafts, and 
Environment, Respectively. The curved solid line is the 
band-pass filtered vertical velocity based on which the 
events were identified. (from Wang 1993). 
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convective updrafts and downdrafts in turbulence transport 


and in the microphysical processes. 


B. GEOMETRIC CHARACTERISTICS OF EVENTS 


The following section analyzes the average sizes of and 
fractional coverage by the updrafts and downdrafts within 
the BL. Results from the conditional sampling method on the 
number of events occurring in each horizontal flight leg and 
the average size of each event will allow determination of 
the geometric characteristics. The mean event size <d> ene 
average horizontal width of an updraft or downdraft) within 
a defined flight leg is calculated as (Wang 1993): 

Yi dua 


<d>u,d= : 9 ; (3) 
Nu, d 





where dis the width of each event, n is the total number of 
events sampled in the leg, and u and d denote updrafts and 
downdrafts, respectively. The event number frequency N (the 
number of events per unit horizontal length) is calculated 
as (Wang 1993): | 

| Nu,d 


Nu, d=— E 4 
T (4) 





where L is the length of the horizontal leg that provided 
the measurements for conditional sampling. The product of 
<d> and N thus defines the proportion of the leg occupied by 
updrafts or downdrafts, and also defines the fractional area 
occupied by the events when the events are randomly 
distributed in space, which is a reasonable assumption for 


the stratocumulus~-topped BL (Nicholls 1989). 
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Figures 6 through 8 show the vertical profiles of the 
mean event size, number Trequency, and fractional area 
occupied by the events plotted as a function of height 
scaled to the BL height (z/zi). The event size is scaled by 
the BL height. The number frequency is presented as Nez, to 
indicate the number of events on a horizontal distance 
equivalent to the BL height. Unfortunately, the scaling by 
Zi has little significance in this study because the BL 
heights of all FIRE flights are similar. 

In general, it is seen that the mean event sizes 
(Figure 6) are between 10 and 30 percent of z;, or about 
200-300 m in the cases of FIRE. In all cases (except for 
Flight 9), the eGnveeeive event size in the upper BL is 
about 0.2 to 0.3 of the BL height, while the near-surface 
event size is between O.20Z2% -and: 0.1523. The vertical 
variation of the event size appears to behave differently in. 
the three well-mixed cases (flights 3, 4, and 6) compared to 
the three clearly decoupled cases (flights 2, 7, and 8). In 
the well-mixed cases, the event size increases with height 
rapidly (below 0.2z:) from the small near-surface value to a 
constant mixed layer value of 0.2z; to 0.32;. In the 
decoupled cases, however, the event sizes hold constant at 
the near-surface value up to 0.4z;. It is likely that the 
entire surface-based layer has a constant value of 0.10z; to 
Os Lozix The results suggest that decoupling indeed 
Suppresses the growth of the convective events. However, 
because of the insignificance of the z, scaling, the fact 
that all events in the surface-based layer are between 0.10 
and 0.15z; does not necessarily mean that the event size in 
the surface-based layer scales well with the BL height in 


the decoupled cases. 
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<dy a>/ Zi <dy g>/Z 
Figure 6. Mean Event Sizes for Updrafts (+) and 
Downdrafts (‘), Scaled by the BL Depth, for the Horizontal 
Turbulence Legs of all Seven Flights. Circled ofr 


uncircled symbols at the same level go together as the 
updraft/downdraft pair for a given leg. The dotted line 
shows the average height of the cloud base. : 
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Figure 6, Continued. 
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The number frequency of events (Figure 7) is smaller in 
the upper BL than in the near-surface layer in all cases. 
‘Compared to the trend of vertical variation of the event 
size, this decrease seems to correlate with ine Ancrease of 
the event size. Therefore, the decrease in number of events 
is a result of the increase in event size. The strongly 
well-mixed case (flight 3) is the only case that does not 
show sizable decrease in number of convective events with 
height. This appears to agree with the findings in Young 
(1988) for the well-mixed clear convective BL. | 

In flights 2, 4, 7, and 8, the downdrafts are larger in 
size and occupy a larger fractional length (Figure 8) than 
the updrafts in the surface-based layer. In chapter II it 
was noted that these flights had large surface buoyancy 
fluxes compared to the in-cloud buoyancy fluxes. It is 
likely that the larger surface buoyancy flux is generating 
concentrated updraft plumes in the near surface layer. 
Randall (1980) and Wang (1993) suggest that the larger size 
of downdrafts implies positive vertical velocity skewness 
(i.e., the updrafts are more vigorous than the downdrafts). 
In flight 8, vertically well-developed cumulus existed 
below, and often extended up into the main stratocumulus 
deck during this flight. This development would tend to 
concentrate relatively strong updrafts into a smaller area 
in the surface layer, explaining the dominance of the area 
by the downdrafts. The opposite is true for the flights 
with small positive (flights 3 and 9) or negative (flight 6) 
surface buoyancy fluxes. Here the updrafts occupy slightly 
more area than the downdrafts. In contrast, because the in- 
cloud buoyancy flux is large for flight 3, the updrafts 
occupy a larger fractional length than the downdrafts near 


the top of the BL. 


34 











Figure 7. Same as in Figure 6, Except for Event 


Frequency. The product N'z; describes the number of events 
in a horizontal length equivalent to the BL depth. 
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Figure 7, Continued. 
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Figure 8. Same as in figure 6, Except for Fractional 


Length Covered by Updrafts and Downdrafts. 
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Figure 8, Continued. 
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The total fractional length covered by convective 
events varies significantly between the different flights 
and between different levels in the BL within each flight. 
The two clearly well-mixed cases, flights 3 and 4, have the 
largest fractional length coverage by convective events, 
between 0.6 and 0.8. Young (1988) reported that the 
fractional area coverage by thermals in the clear convective 
BL is largest near the surface and near the top of the BL, 
and smallest near the middle of the BL. He also states that 
the fractional length coverage is less than 0.5 throughout 
most of the BL, with minimums ranging from 0.36 to 0.43 in 
the mid-BL. Other studies by Manton (1977), Lenschow and 
Stephens (1982), and Greenhut and Khalsa (1987) also report 
smaller values than the ones found here for the FIRE data. 
In flights 2, 3, and 4, and 6 the fractional length occupied 
by the updrafts and downdrafts increases with height through 
the surface-based layer. This trend is also observed by 
Khalsa and Greenhut (1985). These flights also have the 
highest fractional length values within the cloud mixed 
layer. In flights 7, 8, and 9 the fractional length 
decreases with height in the surface layer, and therefore 
the highest fractional length values are in the surface 
layer. Flights 3 -and 7 are the extreme cases in each 
category, where these trends are clearly seen in Figure 8. 
These trends are related to the strength of turbulence at 
each level in the BL. Stronger turbulence means that a 
greater area will be occupied by updrafts and downdrafts. 
For instance, flight 3 has a large convective vertical 
velocity in the cloud layer, implying strong cloud layer 
turbulence, and weak surface layer turbulence (see chapter 
II, Table 2). The opposite structure is observed in flight 
rs 
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IV. THERMODYNAMIC AND MICROPHYSICAL PROPERTIES OF 
CONVECTIVE EVENTS 


A. THERMODYNAMIC PROPERTIES 


Figures 9 through 13 show the vertical profiles of 
several thermodynamic quantities in both the updrafts and 
downdrafts in the horizontal turbulence legs of the FIRE 
flights. The vertical velocity and liquid water are 
averaged over all updraft or downdraft events in the 
turbulence leg, to obtain an event average. The virtual 
potential temperature and water vapor content are band-pass 
filtered to remove the low frequency perturbations, and then 
averaged over all updraft or downdraft events in the 
turbulence leg to obtain an average perturbation from the 
mean. 

The mean vertical velocity in both the updrafts and 
downdrafts for all flights. is fairly symmetrical with height 
in the BL. The vertical variation of the event vertical 
velocity (Figure 9) is consistent with that of the variance 
of vertical velocity (not shown). Here the nocturnal case of 
flight 3 has the strongest convective events as a result of 
strong radiative and evaporative forcing (Wang and Albrecht 
1994). | 

The virtual potential temperature perturbation 
(approximately referred to as buoyancy excess) profiles 
indicate large variability among the flights. Near the 
surface, the updrafts are in general weak but positively 
buoyant while the downdrafts are negatively buoyant. These 
updrafts and downdrafts contribute to the positive TKE 
buoyancy production indicated in Table 2. The exception is 


with flight 6, where the sign of buoyancy perturbation is 
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Figure 9. Mean Vertical Velocity (in m s"“) for Updrafts 
(+) and Downdrafts (*), Scaled by the BL Depth, for the 
Horizontal Turbulence Legs. Circled or uncircled symbols 
at the same level are the updraft/downdraft pair for a 
given leg. The dotted line shows the average height of the 
cloud base. 
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Figure 9, Continued. 
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completely reversed in correspondence with the negative 
surface buoyancy flux (Table 2). Figure 10 shows that the 
cool updrafts can be identified up to 0.35z:. In the upper 
BL, Signatures of the cool downdrafts and warm updrafts are 
clear in cases of sizable in-cloud buoyancy flux, such as in 
flight 3, and recognizable in some flights with weak in- 
‘cloud buoyancy forcing (e.g., flights 2, 6, and 8)... More 
complications are observed here especially near the cloud 
top. In flight 7, however, the downdrafts appear to be 
warmer compared to the updrafts in the cloud layer. It is 
noted in Table 2 that turbulence in this flight is largely 
generated by surface buoyancy flux and wind shear, while the 
in-cloud generation of TKE is very limited. There is thus no 
coherence between the vertical motion and the buoyancy 
excess. | | 

It is known that both radiative cooling and evaporative 
cooling from entrainment mixing may result in cool 
downdrafts. | The presence of evaporative cooling is 
identified by the so-called buoyancy reversal criterion 
(Siems et al. (1990)): 

NO Rhy (3) 
C, 

Here, A denotes the jump across the inversion, L and C, are 
the latent heat of water phase change and specific heat 
capacity of dry air, respectively, 9 is the equivalent 
potential temperature and q is the total water. FOr 
typical BL conditions, k has a value of 0.23. Average 
inversion jump conditions for the FIRE flights calculated in 
McDowell (1999) are used in Equation 5. It is found that (5) 


is satisfied for all seven flights in this study, suggesting 
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Figure 10. Same aS Figure 9, Except 
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Figure 10, Continued. 


46 














that evaporative cooling plays a role in the cool downdrafts 
in the upper BL in Figure 10. | 

In all cases, the difference in buoyancy excess in the 
middie of the BL is very small and may not be significant. 
Buoyancy production in this region is in general the 
weakest, and sometimes may be negative (not shown). 

A somewhat consistent trend is seen in the vertical 
profiles of water vapor perturbation. Figure 11 shows large 
differences between the positive updraft perturbations and 
the negative downdraft perturbations near the surface in 
‘most flights and also near the cloud top in flights 4 and 6. 
We observe smaller differences in the cloud mixed layer as 
the perturbations in updraft and downdraft become smaller in 
magnitude. Near the surface, large positive water vapor 
perturbations are observed in the updrafts near the large 
moisture source of the sea surface, and large negative 
perturbations are observed in the downdrafts as_ the 
descending air parcels become warmer and drier compared to 
the surrounding environment. Minimum values are generally 
seen in the cloud mixed layer away from the cloud 
boundaries. Here, both updrafts and downdrafts are likely 
saturated, so the small difference between the updraft and 
downdraft is due to small temperature differences, with the 
updrafts being slightly warmer. The larger differences near 
the cloud top are likely the result of entrainment, in which 
case the downdrafts will contain drier and warmer air from 
above the inversion. As a result, some of the entrainment 
mixture may become sub-Saturated. Note that in flight 6, the 
updrafts are cool but moist in the stable near-surface 


environment. 
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Figure 11. Same as in Figure 9, Except for Mean Water Vapor 
Perturbation (in g kg”). 
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Figure 11, Continued. 
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Figure 12 shows the vertical profiles of event mean 
liquid water content for the in-cloud turbulence legs. In 
all cases except for flight 2, there is a consistent trend 
of increasing liquid water content with height in both 
updrafts and downdrafts. The difference between the two 
types of dynamical events is also noticeable. Figure 13 
shows the fractional difference between updraft and 
downdraft cloud liquid water content. At most levels in all 
flights the fractional difference between updrafts and 
downdrafts is positive, suggesting more liquid water content 
in turbulent updrafts. The difference is, however, small in 
magnitude, usually less than 5%. Positive differences 
between updraft and downdraft liquid water content appear to 
be the largest near the cloud base in most flights. This 
trend is consistent with the activation and growth of cloud 
droplets in the updrafts at the cloud base. in: 2LAghe. .S; 
the maximum difference is near the cloud top. This is 
likely the reflection of the increased evaporation occurring 
here as a result of the enhanced turbulent mixing and 


entrainment (Wang and Albrecht 1994). 
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Figure 12. Same as Figure 9, Except for Mean Liquid Water 
Content (in g kg”). 
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Figure 12, Continued. 


52 














0.5 
~0.1 -0.05 0 0.05 0.1 








0.5 0.5 , 
0.1 -0.05 0 005 O41 “0.1 -005 0 005 Of 


q Fa 


Figure 13. Fractional difference in Mean Liquid Water 
Content Between Updrafts and Downdrafts (Fg) for the 
Horizontal Turbulence Legs. Fat is calculated as: 
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Figure 13, Continued. 
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MICROPHYSICAL PROPERTIES 


1. Droplet Concentration and Size 


The volume diameter of cloud droplets is calculated as: 


1/3 








where d is the droplet diameter, n is the number of 
droplets, and k is the size bin of the FSSP. Figure 14 
shows the cloud droplet volume diameters for the cloud mixed 
layers of all flights. In general, the droplet sizes 
increase with increasing height from the cloud base. MThis 
is consistent with results from stratocumulus clouds in 
other regions (e.g., Nicholls 1984). Near the cloud top, 
signs of decreased droplet size are observed in some of the 
flights (flights 2, 3, 4, and 6). It is possible that with 
more cloud top data for the other flights the same decrease 
would be seen there as well. 

The largest droplet size is found in flight 6, where 
the maximum droplet diameter near the cloud top is as large 
as 23 pm. This flight also has the lowest BL depth of all 
flights (751 m) and a relatively thick cloud layer depth 
(383 m). The large droplet size is a reflection of the deep 
cloud layer. Droplet diameters of all other flights average 
around 15 um near the cloud top. 

The difference in droplet sizes between the updrafts 


and downdrafts can be clearly identified in Figure 14, where 
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the dot symbols (downdrafts) are in general to the right of 
the corresponding plus symbols (updrafts). The fractional 
difference in volume diameter between the two types of 
convective events is shown in Figure 15. Although small, 
this quantity is found to be consistently negative in most 
cases. This trend indicates that the downdrafts contain 
the larger droplets originating from higher levels in the 
cloud, while the updrafts contain droplets still in the 
growth process. Larger differences are observed near the 
base and top of cloud, with a minimum difference in the mid 
cloud range, which is consistent with droplet formation in 
the updrafts near the cloud base and evaporation due to 
entrainment in the downdrafts near the cloud top. 

The number concentration of droplets in the updrafts is 
consistently higher than in the downdrafts at all levels” 
within the cloud (Figure 16). Since activation of cloud 
droplets is occurring in the updrafts and evaporation of 
droplets is occurring in the downdrafts, these observations 
are reasonable. Nicholls (1984) showed that the number 
concentration of droplets in the stratocumulus cloud is 
fairly constant with height at any given horizontal location 
in the cloud, which is apparent in a sounding through the 
entire depth of the cloud (Figure 3). This cannot be 
observed here since each of the turbulence legs takes place 
at a different horizontal location within the cloud layer, 
which causes the apparent scatter a9) the number 


concentration between the different legs. 
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Figure 14. Same as in Figure 9, Except for Volume Diameter 
of Cloud Droplets (in pm). 
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Figure 14, Continued. 
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Figure 15. Same as Figure 13, Except for Fractional 
Difference for Volume Diameter of Cloud Droplets (Fay). Fay 





is calculated as: (d, ~4,,)/5(d, +d,). 
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Figure 15, Continued. 
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Figure 16. Same as in Figure 9, Except for Number 
Concentration of Cloud Droplets (in cm”). 
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Figure 16, Continued. 


62 











2. Cloud Droplet Spectra 


Many past studies indicate strong vertical and 
horizontal variations of cloud droplet spectra in marine 
stratocumulus due to various reasons. Martin et al. (1995), 
for example, indicated the effects of cumulus penetration on 
cloud microphysics. Noonkester (1984) observed horizontal 
variability in the droplet spectra of a stratus layer off of 
the southern California coast, attributing the variability 
to the horizontal variation in vertical mixing and relative 
humidity. Using large eddy simulations as a tool, Kogan et 
al. (1995) suggest that the cloud microphysical properties 
are significantly affected by BL dynamics, and that the 
cloud microstructure is significantly asymmetric between 
updrafts and downdrafts, especially near the cloud top and 
base. Their results suggest that ‘natural’ horizontal 
Variability, associated with convective events, exists in 
all stratocumulus. However, such natural variability has 
not been analyzed from observations, although it is 
extremely important to the understanding of cloud 
formation/dissipation processes. | 

The integrated microphysical properties, such as the 
liquid water content, droplet volume diameter, and the 
number concentrations, were discussed in the previous 
sections. This section will focus on the variability in the 
droplet spectra. The results will be consolidated through 
consistency among different flights. 

Some consistency is seen in the changes in droplet 
spectra with height within the stratocumulus cloud deck. 
Figure 17 shows an example of the mean droplet spectra in 
the convective updrafts and downdrafts at three 


representative levels in the cloud layer. The results are 
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from flight 6, but are typical of those observed during 
FIRE. The droplet spectra in panels a, b, and c in Figure 17 
were obtained from measurement legs near the cloud base, in 
the middle of the cloud layer, and near the cloud top, 
respectively. 

Figure 17 shows the increase in droplet modal size from 
the cloud base to the cloud top. Integration of the spectra 
will naturally result in the increase in mean droplet 
diameter sHNown in Figure 14. It is seen in Figure 17 that 
the updrafts contain a larger concentration of droplets than 
the downdrafts, with the droplet spectra of the defined 
environment falling in between the updraft and downdraft 
droplet spectra. The difference in peak .droplet 
concentration appears to increase towards the upper cloud 
layer. Near the top, the shape of the droplet spectra in 
updrafts and downdrafts seem to deviate rather 
significantly, with significant reduction in the number of 
small droplets occurring in the downdrafts. The large 
droplets do not appear to be affected... Such change in the 
shape of the droplet spectra result in larger mean droplet 
diameter in the downdrafts, as seen in Figure 14. 

Strong horizontal variability on a scale of tens of 
kilometers is also observed in most of the flights analyzed 
in this study. Figure 18 shows the mean droplet spectra in 
the convective events from a different segment of the same 
cloud top turbulence leg in flight 6.,shown in Figure lic. 
Here, the total number concentration is significantly less 
although the modal droplet size is nearly the same. In 
addition, the downdrafts have a larger droplet number 
concentration than the updrafts, the exact opposite to the 


droplet spectra shown in Figure 17c. The reduction in total 
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Figure 17. Cloud Droplet Spectra from Three in-cloud 


Measurement Legs in Flight 6. The vertical axis 
GN/dlogD, where N is the droplet number concentration 
cm”) and D is the droplet diameter (in wm). a) near 
cloud base, b) near the middle of the cloud, c) near 
cloud top. The symbols for updrafts, downdrafts, and 
defined environment are given in the legend. 
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Figure 18. Same as Figure 17, Except for a Different 
Segment of the Same Turbulence Leg as in Figure lic. 
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droplet number (while retaining the shape of the spectra) 
compared to the other time segment of the same measurement 
leg (Figure lic) appears to indicate inhomogeneous 
entrainment mixing, as discussed in Baker et al. (1980). It 
is not clear, though, whether the change in the size spectra 
in the updraft and downdraft events is typical and what 
caused this. change. With the limitations of in situ 
observations, the answer to these questions may have to come 
from large eddy simulations involving explicit cloud 
microphysics. 

Figure 19 shows the droplet size distribution from two 
different turbulence legs from flight 3 that were flown in- 
cloud near the cloud base. Figure 19a has a spectral shape 
similar to the typical case in Figure 17, yet the plot 
indicates that the droplet concentration is largest in the 
defined environment, slightly higher than the updraft 
concentration. This implies that the turbulence is not as 
active in this region of the cloud base. As a result, the 
turbulence updraft and downdraft may not be well correlated 
with the activation and growth of- the cloud droplets. 
Figure 19b shows a bimodal size distribution, with a larger 
concentration and smaller droplets in the updrafts in both 
modes. This spectral structure indicates that, in this 
well-mixed turbulent BL, there are downdraft droplets being 
recaptured and recycled in the updrafts at this particular 
horizontal location near the cloud base. This phenomenon 
was reproduced in a large eddy simulation model by Kogan et 
ale (1999) < 

We observed horizontal variability in the cloud droplet 
spectra in the cases where cumulus clouds penetrate the 
stratocumulus deck from below. Figure 20 illustrates one 


example of this occurrence in the FIRE flights using 
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Figure 19. Same as in Figure 17, Except for Two 
Horizontal Turbulence Legs Near the Cloud Base in Flight 
3. Legs are between 60 and 85 km apart. 


68 














measurements from flight 7. Figure 20 shows two cloud 
droplet spectra from the same horizontal measurement leg 
near the cloud top. During this in-cloud flight leg, the 
observer’s notes indicate the presence of thiiek 
stratocumulus, with some regions of penetrating cumulus 
observed. The droplet spectra from the first segment 
(Figure 20a) appear similar to those shown earlier for 
flight 6 for the updrafts, downdrafts, and defined 
‘environment, and thus it is consistent with the droplet 
spectra of stratocumulus with no penetrating cumulus (Figure 
17). In Figure 20b, however, the total number concentration 
is much higher, and the updraft and downdraft concentration 
is much higher than that of the defined environment. Also, 
the droplet sizes are smaller than those in Figure 20a. 
This is consistent with the findings of Martin et al. 
(1994), who observed that the intrusion of cumulus clouds 
resulted in a localized increase in droplet concentration 
and liquid water content in the stratocumulus layer, with a 
localized decrease in droplet size as more droplets compete 


for available water vapor. 
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Figure 20. Same as Figure 19, Except for Flight 7 Where 
Penetration of the Stratocumulus by Cumulus was Observed. 
Cumulus penetration was recorded by the in-flight 
observer in (b). 
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on CONVECTIVE MASS FLUX PARAMETERIZATION 


Betts (1975) developed the mass flux parameterization 
originally used to describe vertical flux transport in 
cumulus clouds. Betts (1978) later extended this to 
describe convection in the sub-cloud layer. The mass flux 
parameterization was originally developed to represent the 
vertical flux transport using the properties of the 
convective elements only. The parameterization applied to 
‘Clear convective BL rather successfully (e.g., Khalsa and 
Greenhut (1985) and Greenhut and Khalsa (1987), for the 
trade wind cumulus BL, and Young (1988), for the clear 
convective BL). It 1s not clear, though, whether such 
parameterizations can be applied to the cloudy BL, as the 
convective activity in the cloudy BL is weaker and the 
contrast between the updrafts and downdrafts is in general 
smaller in the cloudy case. 

In this parameterization from Betts (1978), fluxes of 


dry and moist static energy are represented by the product 


of a convective mass flux, @*, and the difference between 


properties inside and outside convective elements, 
fot Wg 
wx =o. (x, —x,) (7) 


where @*, is the convective mass flux associated with the 
quantity x, and xX, and xg are the conditionally averaged 
value of x for updrafts and downdrafts, respectively. Here, 
Equation (7) 1s tested using the observations from the seven’ 


FIRE flights in stratocumulus clouds. The convective mass 


flux, @*, is calculated using the buoyancy flux, water 
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vapor flux, liquid water flux, and ozone flux from the 
horizontal turbulence legs at different levels within the 
BL. The results of the m* calculations are shown in Figure 
Zi The calculations using the ozone flux from flight 9 
were omitted due to a malfunction of the Model Mark II ozone 
sensor during that flight. 

The results show a fair amount of consistency between 
the values of @* calculated from the different fluxes. 
There is some scatter in the data points, the least amount 
of scatter being seen in the well-mixed cases of flights 3, 
4 and 6, where the fluxes throughout the entire BL are 
rather strong. iIn these three flights, there is a trend of 


'increasing @* with height in the lower one third of the BL. 


In the upper two thirds of the BL, @* remains farriy 


constant with the exception of flight 3 near the cloud top. 


Near the surface, @* averages about 0.2 ms, then increases 


to about 0.3 to 0.4 ms’ in the upper BL. This trend for o* 
can also be seen, though not as clearly, in the decoupled 
cases of flights 7 and 8. The results for flights 2 and 9 
have significant scatter at various levels in the BL and do 
not show a clear trend. The values and variation of the 
calculated convective mass flux with height are similar to 
those obtained by and Khalsa and Greenhut (1985) and 
Greenhut and Khalsa (1987), who calculate 0.15<w*<0.25 for 


latent heat, sensible heat, and along-wind momentum fluxes. 
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Figure 21. Convective mass flux, @* (in ms") as a 
Function of Height Scaled by the BL Depth. Fluxes 
corresponding to each symbol are: 9, (+), gq, (*), qi (0), 
and O03 (x). The dotted line shows the average height of 
the cloud base. 
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Figure 21, Continued. 
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V. SUMMARY, CONCLUSIONS, AND DISCUSSION 


A. SUMMARY AND CONCLUSIONS 


‘The -overall. -obyective of (his thesis was to ‘better 
understand the convective nature of the stratocumulus-topped 
boundary layer. This study used data collected by the NCAR 
Electra during seven of the FIRE flights. Of the seven 
flights analyzed in this study, three were determined to be 
well-mixed, three were decoupled, and one was uncertain. 

We used a conditional sampling method to identify the 
turbulent updrafts and downdrafts in the  horazontal 
turbulence legs at different levels within and below cloud. 
We selected minimum event width and vertical velocity 
criteria based on analysis of the data prior to the. 
conditional sampling and also on the criteria used in 
Nicholls (1989), and investigated the geometric 
characteristics and leg-averaged properties of thermodynamic 
and microphysical quantities within the updrafts = and 
downdrafts. 

An in-depth analysis of the geometric properties of the 
convective events revealed some consistent trends. The 
vertical variation in the event size in the well-mixed BL 
cases was observed to increase steadily with height compared 
to the decoupled cases, where event size remained fairly 
constant up to about 0.4z;, suggesting that decoupling 
indeed suppresses the growth of the convective events. The 
number frequency of events was observed to be smallest in 
the upper BL in all flights. This decrease in the upper BL 
seems to correlate with the increase in the event size. The 


trends in the vertical variation of fractional area occupied 


4D 


by updrafts and downdrafts are related to the strength of 
turbulence at each level in the BL. Stronger turbulence 
implies a greater area occupied by updrafts and downdrafts. 
We observed great consistency in the analysis of many 
of the thermodynamic quantities in the updrafts~ and 
downdrafts. The mean vertical velocity in both the updrafts 
and downdrafts for all flights is fairly symmetrical with 
height in the BL Flight 3, the well-mixed, nocturnal case 
containing strong radiative and evaporative forcing, had the 
strongest convective events. Large water vapor perturbations 
were observed near the sea surface (moisture source) and 
cloud top (warm, dry air source), with minimum values 
generally seen in the cloud mixed layer away from the cloud 
boundaries. Both updrafts and downdrafts showed a consistent 
trend of increasing liquid water content with height in the 
cloud layer. The fractional difference in liquid water 
content between updrafts and downdrafts was also 
consistently positive, with the largest positive liquid 
water differences occurring near the cloud base, consistent 
with activation of cloud droplets in the updrafts here. We 
observed less consistency in the virtual potential 
temperature perturbation (or buoyancy excess) profiles, 
which varied greatly among the flights. This variability is 
likely the result of the variation in cloud top and surface 
forcing between the different cases. The average inversion 
jump conditions for the FIRE flights from McDowell (1999) 
were applied to the buoyancy reversal criterion of Siems et 
al. (1990). The buoyancy reversal criterion was satisfied in 
all seven flights, suggesting that the cool downdrafts 
identified in the upper boundary layer are likely a result 


of both radiative cooling and evaporative cooling. 
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Analysis of updraft/downdraft microphysical parameters 
revealed consistent trends as well. In nearly all of the 
flights, the downdrafts contain fewer and larger droplets 
while the updrafts contain more and smaller droplets, 
consistent with droplet formation and growth in the updrafts 
and evaporation due to entrainment drying in the downdrafts. 

Both, vertical and horizontal variability were observed 
in the updraft and downdraft cloud droplet spectra from the 
FIRE flights. As already noted, the updrafts typically 
contain a larger concentration of droplets at all droplet 
size bins than the downdrafts for a given in-cloud leg. 
Additionally, the droplet spectra of the defined environment 
falls in between the updraft and downdraft droplet spectra, 
and there is typically no shift in the peak droplet size 
between the updrafts, downdrafts, and environment at a given 
height. There is a shift to larger droplet sizes with 
increasing height in the cloud, independent of § the 
convective events. The difference in peak droplet 
concentration between updrafts, downdrafts, and environment 
appears to increase near the cloud top. Additionally, the 
number of small droplets occurring in the downdrafts appears 
to decrease near the cloud top, resulting in the observed 
larger mean droplet diameters in the downdrafts here. In 
most of the flights analyzed, we observed significant 
differences in the droplet spectra between two legs at the 
Same vertical level in a particular flight, indicating the 
horizontal variability of the cloud microstructure on a 


scale larger than the boundary layer internal circulation 


(50 to 100 km). The differences seen are likely the result’ 


of localized changes in the extent of turbulent mixing 


occurring near the cloud top and base and also the 


ce 


interaction with cumulus forming beneath the stratocumulus 


deck in a localized area. 


The validity of the mass.flux (@*) parameterization, 
based on methods used in previous work for the clear 
convective BL by Betts (1975 and 1978), Khalsa and Greenhut 
(1985) and Greenhut and Khalsa (1987), was tested for the 
stratocumulus-topped BL. The results suggest that the 
parameterization is applicable to the FIRE data, especially 


in the cases of stronger turbulence where there is less 
scatter among the values of w* calculated using different 


fluxes. Near the surface, @ averages about 0.2 ms", then 


1 in the upper BL, which is 


increases to about 0.3 to 0.4 ms” 
comparable to values obtained by and Khalsa and Greenhut 


(1985) and Greenhut and Khalsa (1987). 


B. DISCUSSION 


This study has looked closely at the convective 
structure of the BL using calculations based on data from 
observations. We do not present error statistics for the 
calculations of the convective event mean profiles, 
geometric properties, and droplet spectra made in this 
study. in one aspect, there is no accepted method of error 
representation Lor these types of calculations. 
Nevertheless, confidence in the results is obtained through 
consistency between results from different flights. | 

The geometric and thermodynamic properties identified 
from this study are also supported by the work of other 
researchers using other datasets. The Current. researcn 
provides, for the first time, the microphysical properties 


in the convective events. These observational results are 
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consistent with those from large eddy simulations with 
explicit cloud microphysics, such as in Kogan et al. (1995). 
The differences in the microphysical properties between the 
turbulent updrafts and downdrafts clearly revealed the 
physical processes occurring in each type of event. More 
well-designed in-cloud measurements are desirable to further 
quantify the findings here. Furthermore, this study clearly 
identified the differences in the event properties between 
the well-mixed cases and the decoupled cases. These results 
should be quantified with more cases in each category to 


achieve high statistical significance. 
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